A NOVEL METHOD FOR MRI BRAIN IMAGE SEGMENTATION 
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khstract This paper presents an automatic segmentation 
algorithm for MRI (magnetic resonance imagind brain image. The 
algorithm integrates morphological Algorithm and histogram- 
based fuzzy Cmeans(FCM) to segment brain into three major 
classes of gray matter (GM), white matter (WM) and 
cerebrospinal fluid (CSF). First, morphological is used to Removal 
of extracranial tissues to get brain image and then histogram- 
based FCM algorithm is used to segment the imagesThe 
algorithm was implemented in MATLAB code.Satisfied results 
were obtained. 
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I. INTRODUCTION 

When working with medical images it is oftea of interest tD 
delineatB interesting areas or vouimes.The prc)cess of finding 
those is called segmeatationMagnetic resonance images (MRI) 
of human biBin typically contain three tissue classesigray 
matter (GM), v\*ite matter (WM) and cerebrospinalfluid (CSF), 
as shown in Figure ( 1) . 

siiuiple MR! : CSF TPM: Grey-matter TPMr White-matter TPM: 



Fig.l three tissue classes:gray matter (G M ), white matter (W M ) and 
cerebrospinalfluid (CSF) 

Fully automatic, accuiate and robust biBin tissue 
classi" cation from these MRI data is of great importance for 
research and clinical studies of the normal and diseased human 
brain segmentation methods for performing tissue 
classification are mainly hindered by the following three 
imaging artifacts. First, because the intensities of MR images 
are sensitive to the conditions or imaging parameters of the 
scannei:s, we don't have a dear-cut correspondences between 
the intensities and the biological tissues. There may have big 
oved^DS of the intensity distribution ranges among different 
tissues. Second, the second difficult is the inhomogendties of 


the magnetic fidd and biological variations in different 
stmctures bdonging to the same tissue class. Third, the last 
main problem is the partial volume effects in MR images 
which made the tissue boundary blurred[l]. Because of these, 
using conventional intensity-based method to classify MR 
images is problematic. In the p^Der, we proposed an algorithm 
not only taking use of the characteristic of the intensity 
distribution but also the morphological character 

II TheMORPHOLOGICAL METHOD 

1) Background on mathematical morphology: Image 
morphologyprovides a way to incorporate neighborhood and 
distance information into algorithms (see Serra, 1982; 
Haralick et al., 1987 for detailed treatment of morphological 
operators). The basic idea in mathematical morphology is to 
convolve an image with a given mask (known as the 
structuring element) and to binarize the result of the 
convolution using a given function.Choice of convolution 
mask and binarization function depend on the particular 
morphological operator being used.Binary morphology has 
been used in several segmentation Systems. 

• Erosion: an erosion operation on an image 1 containing 
labels and 1, with a structuring element S, changes 
the value of pixel i in 1 from 1 to 0, if the result of 
convolving S with 1 , centered at i , is less than some 
predetermined value. We have set this value to be the 
area of S, which is basically the number of pixels that 
are 1 in the structuring element itself. The structuring 
element (also known as the erosion kernel) determines 
the details of how a particular erosion thins boundaries. 

• Dilation: dual to erosion, a dilation operation on an 
image 1 containing labels and 1, with a structuring 
element S, changes the value of pixel i in 1 from to 1, 

if the result of convolving S witii 1 , centered at i , is more 
than some predetermined value. We have set this value to be 
zero. The structuring element (also known as the dilation 
kernel) determines the details of how a particular dilation 
grows boundaries in an image. 

• Conditional dilation: a conditional dilation is a dilation 
operation with the added condition that only pixels that are 1 
in a second binary image, Ic, (the image on which the dilation 
is conditioned), will be changed to 1 by the dilation process. It 
is equivalent to masking the results of the dilation by the 
image Ic. 
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• Opening: an opening operation consists of an erosion 
followed by a dilation with the same structuring element. 

• Closing: a closing operation consists of a dilation 
followed by an erosion with the same structuring element.[2] 

2) Removal of extracranial tissues to gethrain image: 
An important pre-processing step is the removal of 
extracranial tissues because our segmentation technique labels 
the entire brain into three classes: gray matter (GM), white 
matter(WM) and cerebrospinal(CFS). In the paper, 
morphological method were used to remove extracranial tissue. 
The results for one MRI image shown in Figure. 2. 



Fig.2 . (a) Original image 


(b) after morphology 


lII.The Fuzzy C-means Method 

use of fuzzy segmentations over hard segmentations in 
defining anatDnlcal structuies. Fuzzy segmentations retain 
more information from the original image than hard 
segmentations by taking into account the possibility that moro 
than one tissue class may be present in a single voxel and has 
previously been used with some success in the fuzzy 
segmentation of magnetic resonance (MR) as well as for the 
estimation of partial volume [3].It clusters data by computing a 
measure of membership ,called the fuzzy membership, at each 
voxel for a specified number of classes. The fuzzy membership 
function, constrained to between zero and one, reflects the 
degree of similarity between the data value at that location and 
prototypical data value or centroid, of its class. Thus, a high 
membQ:ship value near unity signifies that data value at that 
location is "close " to the centroid for that particular class. 

Ij Histogram-hased FCM a/gon't/im.-Histogram-based FCM 
method was used in the paper for segment the brain, which only 
use the I-D histogram and thus can greatly simplify the 
segmentation computation. The algorithm based on the 
principle that for a gray level only belongs to two neighbor 
parts with different membership, and the nearer with the 
centroid of that part, the higher membership it has, thus greatly 
simplify the fuzzy segmentation computation. We get good 
segmentation results by optimizing the objective function of 
fuzzy Restrained FCM clustering. Although FCM algorithm 
isrobust to initial centroid, proper selection will generally 
improve accuracy and convergence of the algorithm. In order to 
closely and easily estimate the initial centroids,we developed a 
weighted histogram with very similar peak and vale property 
and smoother than the original histogram. 


2) Initial Centroids: In order to properly select the initial 
centroid values, a weighted histogram PH (k) (including the 
background ) was constructed by using gaussian 

function as smoothed kernel estimator of image histogram 
(see equation 1) and used to estimate the three peaks as the 
initial centroid values estimation vj , k=l,2,3, for 

WM,GMandCSF. 


1 255 

Pmax =max(P^(5)) 


-ig-kf 


(1) 


3) The procedure to find the three peaks is as follomngs: 
1 : find the gray level corresponding the maximum peak of 
weighted histogram. 

2. find the second large peak: a simple trick that frequently 
works well enough is to look for the second peak by 
multiplying the histogram value by the square of the distance 
from the first peak. This gives the preference to peaks that are 
not close to the maximum. So, if tiie largest peak is at level j 
in the weighted histogram, select the second peak asl4]: 

msx{({k-j)'P^{k)p<k<255)} ^2) 

In the same may, we can find the gray level corresponding 
to third largest peak. Sort order the three gray levels by 
ascending, we get the corresponding centroid values.Figure.3 
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Fig.3.(a) Original grayscale (b)Weighted histogram 

image histogram 
Develop the Restrained HistogrBm Fuzzy Men±)ership 
Function: 
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K : the gray level of the image 
i : the classes to be segmented 


(3) 
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Compute the objective function JHFCM of histogram based 
FCM method. 

255 3 ^ 

j™=i:E".>.^-v||H(k) 

k=0 z=l (4) 

By minimizing the objective function, we get the 
Optimized ^^ and ' . 
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13x256 


fuzzy histogram clustering division matrix 


x3 the ceatroids matrix of difeteit tissues in brain 
Properiy selecting weighting exponent m is also important 
for producing good cluster in FCM algorithm. Figure.4 
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Fig.4.a)Restrained histogram fuzzy membership function 


b) fuzzy membership function of image 

Though there is no theoretic methodology, experiment 
results deoDnstrate its ffectiveQess and arrive at the conclusion 
that the optimal selected range m is [ 1 .5 2.5].In the p^Der, we 
usem=1.75 

■(")_V("-i)||<e 

II n.. iterative times 

,fhen using potential function as height readjustment 
function to modulate the clustering division matrix and get 
optimized thresholds for segmenting image based on maximum 
cluster rule. 


F,.(k) = p^(vJ[/,.(k) 

T. =x. , F.(x.) = F.^i(x.) i=l2 


(6) 


we use this method to segment brain( image(a) in figure 
2)and get the threshold Tl = 25(for CSF , CSF and the 
background are in the same class but our algorithm segmented 
CSF from the background , therefor it is no necessary to 
classify them) , T2= 112 and T3=146. results is shown in 
Figure. 5. 



Fig.5.Results of our segmentation 

IV. Conclusion 

This paper described an MRl brain image segmentation 
algorithm which works in two steps: morphological method 
and histogram-based fuzzy C-means (FCM) for segmenting 
the brain image into three tissues GM,WM and CSF. Satisfied 
results were obtained. 
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